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Abstract 

Using our recent proposal for denning gauge invariant averages we give a general- 
covariant formulation of the so-called cosmological "backreaction" . Our effective covari- 
ant equations allow to describe in an explicitly gauge invariant form the way classical 
or quantum inhomogencitics affect the average evolution of our Universe. 



1 Introduction 



It is well known that the homogeneous and isotropic Friedman-Lemaitre- Robertson- Walker 
(FLRW) metric describes the geometric properties of our Universe only on sufficiently large 
scales of distance, and has to be interpreted as the "averaged" cosmological metric: namely, 
the metric emerging from an appropriate smoothing-out of the local inhomogeneities and 
anisotropics. The same interpretation of averaged variables has to be assigned to the matter 
energy and momentum density, sourcing the FLRW metric in the cosmological Einstein 
equations. 

A problem (a rather old one, see for instance [1]) thus appears due to the fact that 
the Einstein equations for the averaged geometry are different, in general, from the aver- 
aged Einstein equations. This is because the averaging procedure does not commute, in 
general, with the non-linear differential operators appearing in Einstein's equations. As a 
consequence, the dynamics of the averaged geometry is affected by so-called "backreaction" 
terms, originating from the contribution of the inhomogeneities present in the metric and 
matter sectors. 

Interest in these themes has considerably risen after the suggestion that the cosmic 
acceleration, recently observed on large scales, could be unrelated to phantomatic dark- 
energy sources, but - perhaps more simply - to the dynamical effects of the backreaction 
(see e.g. [2]) thus solving the well-known "coincidence problem" . This possibility has focused 
current researches on a new problem: the gauge invariance of the averaging procedure. 
Indeed, in the absence of gauge invariance, the computed backreaction effects depend not 
only on the hypersurface chosen to compute the average integrals, but also on the chosen 
coordinate frame (see e.g. [3]- [6] for recent discussions). 

At present, the most commonly used averaging procedure in a cosmological context is 
based on the foliation of space-time into three-dimensional hypersurfaces comoving with 
the matter sources, and on the volume integration over such spacelike hypersurfaces [7] 
(see [8] for a recent review, and [6, 9] for an extension to more general hypersurfaces). 
This procedure is applied, in particular, to the scalar part of the cosmological Einstein 
equations, i.e. to the so-called Hamiltonian constraint and Raychaudhuri's equation. By 
spatially averaging such equations on a comoving domain D (and considering, for simplicity, 
dust fluid sources) one obtains [7]: 

(^f-'-f^-lmn + {K M , (1.1) 

-£ = *£v)D-\m. (i.2) 

an o o 
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Here the dot denotes derivatives with respect to the cosmic time of the comoving syn- 
chronous coordinates, op is an effective scale factor, related to the volume Vd of the 
integration domain (normalized by a reference volume scale Vd ) by an = (Vb/Vb,,) 1 / 3 , 
where 

V D = J ^xy/\det gij \, (1.3) 

and gij is the intrinsic metric of the comoving hypersurfaces. Also, the brackets denote 
spatial average over D, namely 

(• ■■)d = ^ j D d Z x^j\&et 9ij \ (•••), (1.4) 

p is the energy density of the dust sources, and 1Z is the scalar intrinsic curvature associated 
with the spatial metric gij. Finally, (Q) is a correction called "kinematical" backreaction, 
arising in a particular gauge (see section 2) from the averages of two scalar quantities: the 
trace of the expansion tensor and the scalar shear a 2 : 

(Q)D = l((0 2 ) D -(e)l)-2(a 2 ) D . (1.5) 

The above averaged equations are obtained within a spatial slicing of the space-time 
manifold induced by the flow lines of the matter sources. In this paper we present a 
covariant version of the effective equations for the averaged cosmological quantities (called 
hereafter, for simplicity, averaged cosmological equations), based on a recently proposed 
gauge invariant averaging prescription [4] : the backreaction effects we obtain depend on the 
hypersurface chosen to define the physical observer, but do not depend on the particular 
choice of coordinates. In the appropriate class of gauges we recover the recent results given 
in [6, 9]. In addition, for an observer at rest with respect to the matter sources, we recover 
the same results as in [7, 8] . 

The paper is organized as follows. In Sect. 2 we present an explicitly covariant version 
of the Hamiltonian constraint, Raychaudhuri's equation, and of the projected conservation 
equation for a generic foliation of spacetime and energy-momentum tensor. In Sect. 3 we 
briefly summarize the main aspects of our gauge-invariant averaging prescription, and we 
derive the corresponding general-covariant version of the averaged cosmological equations. 
Our conclusive remarks are briefly presented in Sect. 4. 

2 Covariant ADM equations 

In order to introduce spatial averages of physical quantities we consider a general class of 
foliations of spacetime by spacelike hypersurfaces T,(A) over which a scalar field A(x) takes 
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constant values (the so-called level-sets of A). Let n p be the future-directed unit normal to 
E(A), defined by 

= ' " x = _1 (2 - 1) 

(we are using the metric signature (—,+,+,+)). Let us also introduce the projector h pu 
into the hypersurfaces by: 

hfiv = g^u + n^riu , h w h p v = h^ u , h^rf = . (2.2) 

The Einstein equations G pv = T pv (we use units in which 8ttG = 1) can then be projected 
along and h? v , and give rise to three (sets of) equations that can be chosen in the 
following form: 

G^n v = T pv n»n v = e, (2.3) 
G^h; = T pu n^ p = J p , (2.4) 

= T^tfh" - h pa T = S pa - h pa T . (2.5) 

They correspond to an explicitly covariant version of the so-called Arnowitt-Deser-Misner 
(ADM) equations. 

It is always possible to make contact with the more conventional ADM formalism by 
choosing a class of gauges in which the scalar field A(x) is homogeneous. We shall call it 
the ADM gauges. In such a gauge the normal vector n p takes the form: 

n M = iV(-l, 0,0,0) , n M = l(i,_AP), (2.6) 

where N and N l are, respectively, the so-called lapse function and shift vector. In this 
gauge: 

K = Sf, h° = 0, h ij =g ij , (2.7) 

where g^ is the induced 3-metric (or first fundamental form) on S. The spacetime metric 
of the foliated spacetime is given by: 

ds 2 = g pu dx^dx u = -N 2 dt 2 + g ij (dx i + N i dt)(dx j + N j dt) , (2.8) 

and its inverse by: 

□ = g ^V p V u = -N~ 2 (V - N*Vi) 2 + Wg ij ViVj . (2.9) 

where ^g % i is the inverse of the 3x3 metric gij . 

In this class of ADM gauges Eqs. (2.3) and (2.4) reduce to the convential form of 
the Hamiltonian and momentum constraints, respectively, while Eq. (2.5) generates the 
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second order evolution equations. In the ADM context such evolution equations are splitted 
into twice as many equations, those defining the second fundamental form (or extrinsic 
curvature) Kf, and those describing the evolution of K\ itself (see, e.g. [8]). 

In order to give a covariant and gauge invariant formulation of the cosmological back- 
reaction we shall make use of the covariant Eqs. (2.3)-(2.5), which lead to an explicitly 
scalar form of the Hamiltonian constraint and Rachayduri's equation. To this purpose, let 
us first consider the spacetime flow generated by the timelike vector field n M , and define the 
projected expansion tensor of the flow worldlines as 

= h^\7 a n fj = -h^Q + o> + , (2.10) 

where 

e = V M n^, op, = h a ^ u (v (a n p) - hi a pV T n T ^j , = h^V [a n p] , (2.11) 

are the expansion scalar, the shear tensor and the rotation tensor, respectively. In our case, 
with n M given by Eq. (2.1), we have a zero rotation tensor and we can write 

where a 2 is the shear scalar. 

Consider first the Hamiltonian constraint (2.3), whose left hand side can be decomposed 
as follows: 

R^n" + \r = \ (e 2 - e^) + l -n s . (2.13) 

In the ADM gauge the term in brackets on the right hand side can be expressed in terms 
of the extrinsic curvature as: 

e 2 - e£6£ = K 2 - K)K{ , (2.14) 

while 1Z S is a scalar which, in the ADM gauge, goes over to the intrinsic scalar curvature 
1Z associated with the induced metric gij. We can thus rewrite the Hamiltonian constraint 
completely in terms of scalar quantities as: 

n s + e 2 - e^e^ = n s + ^e 2 - 2a 2 = 2T^n u = 2e . (2.15) 

Coming now to the explicitly scalar form of Rachayduri's equation, we note that it 
corresponds to the linear combination of Eq. (2.3) and of the trace of Eq. (2.5), leading to: 

R^n" = T^hT - \t. (2.16) 
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After some straightforward calculation the above equation takes the form: 

-n^e = 2a 2 + ^e 2 - V>^V M rv) + (T^ - n»n u 

= 2a 2 + ^Q 2 - VKV^j + e + ^T. (2.17) 

So far we have made no assumptions on the form of T^. Let us now restrict ourselves 
to the case of a perfect fluid with: 

T^ v = (p + P^Uv+pg^, (2.18) 

where is the 4-velocity comoving with the fluid, and p and p are, respectively, the 
(scalar) energy density and pressure in the fluid's rest frame (the generalization to several 
non- interacting fluids is straightforward). We stress that and are in general distinct: 
the former depends on the properties of the matter sources, the latter depends on the 
choice of the hypersurfaces on which we want to average, hence should be determined by 
the particular problem at hand. With this model of sources we can then express the basic 
quantities entering our equations (2.15) and (2.17) as: 

e = T^n» = (p + p)(u^) 2 -p, T = T» = -p + 3p. (2.19) 

For later use, we note that the trace T of the energy momentum tensor is unaffected 
by the possible "tilt" (misalignment) between and n M , while the sources of the covariant 
ADM Eqs. (2.3)-(2.5) - namely the objects that we may call the ADM energy density e, 
the ADM pressure tt = 5^/3, and the ADM current J v - become: 

£ = p _ (/0 + |) ) (1.(^)2) j (2 20) 

tt = p-I( /9 + p )(l-(n^) 2 ) , (2.21) 
Jv = {p+p)u p np{l+u^n^)u u . (2.22) 

On the other hand, a straightforward calculation in the ADM gauge leads to: 

(u^f = 1 + < > 1 , (2.23) 

meaning that this quantity is always larger than 1. We can thus introduce a "tilt angle" 
such that sinh 2 ax = (u^n^) 2 — 1, and we can rewrite e and tt in the more convenient form 

e = p + (p + p) sinh 2 ar , tt = p + -(p + p) sinh 2 a-r . (2.24) 

Let us finally consider the condition following from the projected energy-momentum 
conservation law, n^V^T^ = 0. It is known that such condition is not implied by the 
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two projected scalar Einstein equations considered above, and has to be added to the set 
of averaged cosmological equations as an additional independent constraint [7, 8]. In the 
general context we are considering (with / u^), the projected conservation equation can 
be written explicitly as follows: 

[(P + P) u p n p \ + n^p +(p + p) [V^u p n p - e"%u„] = 0. (2.25) 

3 Covariant averaged equations 

Let us begin from the four dimensional integral of a scalar S(x) as defined in [4]: 

I{S,(l)= f d 4 xJ-g(x) S(x) = f d 4 xJ-g(x) S(x)W n (x). (3.1) 

The integration region C .M4 is defined in terms of a suitable scalar window function 
Wn, selecting a region with temporal boundaries determined by the space- like hypersurfaces 
T,(A) (defined in Sect. 2), and with spatial boundary determined by the coordinate condition 
B < ro, where B is a (positive) function of the coordinates with space-like gradient d p B, 
and tq is a positive constant. As we are interested in the variation of the volume averages 
along the flow lines normal to we choose in particular the following window function: 

W n (x) = n"V^(A(x) - A o )0(r o - B(x)), (3.2) 

where 6 is the Heaviside step function, and n M is defined in Eq. (2.1). 

As discussed in [4], if B{x) is a scalar function the integral (3.1) is not only a scalar 
under general coordinate transformations but is also gauge invariant (to all orders): namely, 
it is invariant under the local field reparametrizations induced by any coordinate change 
when old and new fields are evaluated at the same space-time position. If B is not a scalar, 
the spatial boundary can be a source of breaking of covariance and gauge invariance. In [10] 
we will discuss in more detail such a breaking, and confirm that it goes away in the limit of 
large spatial volumes (with respect to the typical scale of inhomogeneities) [4]. Using the 
window function (3.2), the integral (3.1) becomes: 

I(S,A ) = f d A xJ-g(x)5(A(x) - A )(-d^Ad p Af/ 2 e{r Q - B(x)) S(x) . (3.3) 

Let us now consider the derivative of I(S, Aq) with respect to Aq, a quantity that, like 
/ itself, is covariant and gauge invariant (apart from a possible gauge dependence induced 
by the spatial boundary): 

dI{ Q A ^ ] =~J d 4 xj^)6' (A(x) - A ) {-d,A d^Afl 2 e{r G - B{x)) S(x) 
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= -J d 4 x^/-g{x)d 5(A(x) - Aq) [d^x)]- 1 (-d^Ad^A) 1 /^^ - B(x)) S{x). 

(3.4) 

Within adapted ADM coordinates, where A is homogeneous, we can always choose those 
with vanishing shift (i.e. with goo = — N 2 and goi = 0). In such coordinates we have: 

dI( g A y =~J d 4 xV=g~ d 5 (A(t) - A ) (- g M)V*e(r - B{x)) S(x) 

= J d 4 x5(A(t)-A )d [V^(-g 00 ) 1/2 0(ro-B(x))S(x) 

= J d 4 x yJfr\6(A(t) - Aq) [e(r - B(x)) (NQ S + d S) - S(r - B(x))Sd B , 

(3.5) 

where 7 = det^j, and we have used that, in these coordinates, 

e = iV- 1 dologV7- (3.6) 

The above equation can be easily recast into the following covariant and gauge invariant 
form, 

(3.7) 

that reduces to Eq. (3.5) in the special coordinates we have been using. Note that the last 
factor of the above equation is absent if n M is orthogonal to the gradient of B (n^d^B = 0), 
namely if B does not depend on the time coordinate of the gauge used in Eq. (3.5). We 
shall restrict ourselves to this case hereafter. 

Let us define now the covariant averaging prescription for a scalar S(x) on the hyper- 
surfaces of constant A, following [4], as: 

{S)A ° = J^Ao)- (3 - 8) 
Using Eq. (3.5) we can easily obtain the derivative of the averaged scalar (S)a as: 

and Eq. (3.7) immediately gives us the generally covariant version of the above equation: 

d{s) A() _ / d,Ad»s \ I se \ _ I 9 \ 

dA \d,Ad»A/ A J\(-d»Ad»Ay/i/ Ao W^Xi-dpA&Ay/*/^ ^ Uj 

This is the covariant and gauge invariant generalization of the Buchert-Ehlers commutation 
rule [11], and will be the starting point for our generalization of the averaged cosmological 
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equations. However, let us first illustrate the precise connection to the special version of 
this rule obtained in [11]. 

In ADM coordinates, with A homogeneous, Eq. (3.3) reads: 

I(S,A )= [ d 3 x^/\ 7 (t ,x)\ S(t ,x)6 (r - B(t Q , x)) , (3.11) 

where we have called to the time when A(t) takes the constant values Aq, and the averages 
are referred to a section of the three-dimensional hypersurface where A(x) = Aq. This 
is exactly the type of spatial integrals used in [7] (with a domain D determined by the 
condition B(x) < ro), and thus leads to the same definition of averages (see Eq. (1.4)). Let 
us compute, in this case, the corresponding version of Eq. (3.10). Multiplying both sides 
by dtAo we obtain the equation 



d t (S) Ao = [d t S-N*d i S) A +{SNe) Ao -(S) Ao (N@) Ao 

= (Nn^S) Ao + (SNQ) Ao - (S) Ao (NQ) Aq , (3.12) 

which generalises the commutation rule given in [11] to the case of non- vanishing shift 
vector. For N l = the standard result is recovered. 

Let us come now to the covariant formulation of the averaged cosmological equations. 
Starting from the generally covariant volume integral 

1(1, A ) = J d 4 xV^5(A(x) - A )^-d ll Ad^A9(r - B{x)) (3.13) 

we define, along the lines of the previous approach [7], an effective scale factor a such that 

1 da 1 0/(1, A ) 



adA 3/(1, A ) dA 
We then find, using Eqs. (3.7), (3.8), 

1 da _ 1 / 6 
U~dA Q ~ 3 \ (-d^Ad^A) 1 / 2 



(3.14) 



(3.15) 



We are now in the position of presenting the covariant generalization of the averaged equa- 
tion (1.1). Taking the square of the previous equation, using the Hamiltonian constraint in 
the form of Eq. (2.15), and explicitly reintroducing Newton's constant in the formulae, we 
easily obtain: 

T^\ 2 8ttG / e \ 1 / TZ S \ 
.ddA~J ~ ~T~ \ / Ao ~ 6 \ (—d^Ad^A) / M 



1 / a 2 \ 

+ 7T 
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In order to arrive at the covariant generalization of the second Buchert's equation (1.2) 
we start with the simple relation 

1 d 2 a d /l da \ f \ da \ 2 

Using Eq. (3.15), and the general commutation rule (3.10) for the first term on the right 
hand side, we obtain 

Id 2 a 2 I 6 \ 2 1 / d^Ad^(d u Ad v A) \ 

adA 2 ~ 9 Y-dpAdvA) 1 / 2 / + 6 \ (-d^AWA) 5 / 2 " ) 

1/ ® 2 \ l/_d t Ad^_\ 

3 \ (-9,^4) / Ao + 3 \ (-dpAWA)*/* / Ao ■ 1 • J 

Inserting then the covariant Raychaudhuri's equation (2.17) in the last term of the above 
equation we are lead to the final result, 



1 d 2 a 4ttG / e + 3tt \ 1 / V u (n»V ^n u ) \ 1 I df,Ad^{d v Ad v A) 

ZdA 2 ~ "l"\(-^M)/ Ao "3\(-^M)/ Ao + 6\ (-dpAfrA)*/* 

2 



e 2 \ / e 



2 / a 2 \ 



{—dfxAd^A) ) Ao \ {-d.Ad^Af/ 2 /J + 3 \ (-d,Ad»A) j ^ ' (3 - 19) 

where we have used the relation p — 3p = e — 3tt (see Eq. (2.24)). Equations (3.16) and 
(3.19), together with the averaged conservation equation discussed below, are the main 
results of this paper. 

Let us now observe that Eq. (3.16), written in the ADM gauge, and multiplied by 
(d t A ) 2 , reduces to the equation recently presented in [6, 9]: 

(N 2 p) Ao + — (N 2 (p + p)smh 2 a T ) a - -{N 2 K)a 



adtj 3 v 110 3 \ vr /Ao 6 

" l((N 2 & 2 )A Q -(Ne) 2 A) )+l(N 2 a 2 ) Ao . (3.20) 

We have used Eq. (2.24) to replace the ADM parameter e with the fluid proper energy 
and pressure. When sinh 2 ar = (namely when the averaging hypersurfaces coincide with 
those orthogonal to the fluid velocity, = u^), we exactly recover the corresponding 
Buchert's equation (see the second paper of [7]). In the synchronous gauge (N = 1) and 
for dust sources (p = 0) we recover instead Buchert's Eq. (1.1), apart from the additional 
contribution arising from a nonvaninshing tilt angle olt- 

Consider now the second equation (3.19), and let us first note that 

ld 2 a fdAv\ 2 ld 2 a 1 d 2 A da 



dt 2 \ 8t J ad A 2 a dt 2 0A ' 
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We go then to the ADM coordinates, as before, imposing in addition the convenient gauge 
choice N l = 0, and we apply Eq. (2.24) to express e and tt in terms of p and p. By using 
Eq. (3.19) to eliminate the first term on the right hand side of Eq. (3.21) we obtain 

- I {(n 2 Q 2 ) Aq - (ng) Ao ) + \{N^) Ao 

in agreement with [6]. Again, for sinh 2 olt = 0, we also recover the corresponding Buchert's 
equation (see for instance [8]). We may note that, when sinh 2 ay 7^ 0, the "tilt effects" give 
a negative contribution (assuming p + p > 0) to the average cosmic acceleration. 

Let us conclude this section by noting that, as anticipated in Sect. 2, the set of aver- 
aged cosmological equations has to be complemented by the general-covariant average of 
the conservation equation (2.25). Such an operation can be performed straightforwardly, 
according to the general procedure outlined above. We shall present here, for simplicity, the 
covariant version of the averaged conservation equation in the particular case in which the 
space-time foliation is referred to the fluid comoving frame. In such a case, setting n M = 
in Eq. (2.25), and applying the general-covariant commutation rule (3.10), we obtain 

^ {pUo = " {(-d.AdlA)^) Aq ~ {p) ({-dJ»>A)w) ^ (3 - 23) 

(note that, in this case, the scalar A corresponds to the velocity potential of the fluid 
sources). Going, as before, to the ADM coordinates, and multiplying by dtAo, we recover 
the known result 

^{p) D + (NQp) D - (p) (NQ) D = 0, (3.24) 
already presented in various papers [7, 8]. 



4 Conclusion 

The main results of this paper are the covariant and gauge invariant formulation of the 
Buchert-Ehlers commutation rule, Eq. (3.10), and of the effective equations for the averaged 
evolution of a perfect fluid-dominated Universe, Eqs. (3.16), (3. 19), (3. 23). The average is 
performed over a generic class of hyper surf aces, not necessarily orthogonal to the fluid flow 
lines. We stress that our results allow to compute averaged quantities in a completely 
arbitrary coordinate system. 
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The results obtained in this paper can be directly applied to the case of the quantum cos- 
mological backreaction, by using the correspondence between quantum expectation values 
and classical averages performed over all three-dimensional space, as illustrated in details 
in [4]. Hence, in particular, can be applied to study the effect of the quantum fluctuations 
within the canonical formalism of cosmological perturbation theory. 

At the same time, the classical averaged equations (possibly extended to light-like hy- 
persurfaces) may provide a covariant starting point for determining whether present inho- 
mogeneities can significantly contribute to the observed cosmic acceleration. 
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